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Abstract 

In this paper we study, both with theoretical and experimental approach, the effect of iron 
doping in zirconia. Combining density functional theory (DFT) simulations with the experimental 
characterization of thin films, we show that iron is in the Fe 3+ oxidation state and accordingly 
that the films are rich in oxygen vacancies (V"). V" favor the formation of the tetragonal phase 
in doped zirconia {ZrOi'.Fe) and affect the density of state at the Fermi level as well as the local 
magnetization of Fe atoms. We also show that the Fe(2p) and Fe(3p) energy levels can be used as 
a marker for the presence of vacancies in the doped system. In particular the computed position of 
the Fe(3p) peak is strongly sensitive to the V" to Fe atoms ratio. A comparison of the theoretical 
and experimental Fe(3p) peak position suggests that in our films this ratio is close to 0.5. 

Besides the interest in the material by itself, Zr02'-Fe constitutes a test case for the application 
of DFT on transition metals embedded in oxides. In ZrOi'-Fe the inclusion of the Hubbard U 
correction significantly changes the electronic properties of the system. However the inclusion of 
this correction, at least for the value U = 3.3 eV chosen in the present work, worsen the agreement 
with the measured photo-emission valence band spectra. 

PACS numbers: 68.55.Ln,71.15.Mb,75.50.Pp 
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INTRODUCTION 



In dilute magnetic semiconductors (DMS) magnetic impurities, usually transition metals 
(TM), are introduced to produce a magnetic ground state. These systems have been exten- 
sively investigated since the discovery of carrier induced ferro-magnetism in (In, Mn)As^- 
and (Ga, Mn)Aa^, and are believed to be fundamental to fabricate spin-based electronic 
devices. The understanding of DMS physical properties constitutes a challenge for the 
theory as the fundamental mechanism leading to ferromagnetic interaction can be hardly 
explained^. Also experimentally the inclusion and the influence of TM doping is not clearly 
understood. Indeed, while several DMS were predicted to have a Curie temperature (T c ) 
above room temperature, no experimental report of T c > 300-fT has been left unchallenged 
by other studies^. Moreover some results suggest that magnetic impurities, at least at very 
low doping concentration, act as paramagnetic centers 5 . Recently a new class of DMS, based 
on oxides such as zirconia (ZrO-i) and hafnia (HfO-z), has received great attention, after 
the experimental reports of room temperature magnetism in Fe doped HfO^ and ZrO^— 
and the theoretical prediction of high T c in TM doped ZrO-r^-. 

For a better understanding of the magnetic properties of the system, a clear picture 
of its structural and electronic properties is fundamental. As opposite to standard bulk 
materials, where usually the most stable configuration can be unequivocally identified, in 
DMS the inclusion of the dopant can induce stress, disorder and defects in the system with 
many possible configurations close in energy. From one side, theoretically, the modeling 
of the material, also at the first-principles levels, requires some assumptions on the initial 
structure and on the position occupied by the dopant. From the other side, experimentally, 
stress, disorder and defects make difficult to provide a unique interpretation to the features 
observed. Thus a combined approach is the best option. 

Among the structural defects of dilute magnetic oxides (DMO), oxygen vacancies (V") 
are believed to affect the magnetism 8,9,14 . Indeed it has been suggested that delocalized 
electrons, associated with V" , can play a crucial role in the magnetization mechanisms of 
DMO 3 . However, in this model, V" are assumed to always induce delocalized states, which 
can mediate the magnetic interaction. This assumption is true in the undoped oxide, while 
in presence of doping should be verified case by case. 

In the present paper we describe the structural and electronic properties of iron doped 
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zirconia [ZrOi'-Fe] focusing our attention on the role of V" and on their relation with the 
dopant. The role of V" in ZrO^-Fe, and more in general of ZrO% doped with valence +3 
elements {X +3 , with X = Fe,Y, etc..) has been, in part, explored in view of different 
applications, for oxygen sensing^— and more recently for resistive switching memories^ 1 ^. 
For ZrOi'-Fe in particular only few experimental reports exist. Also for TM doped oxides in 
general, no systematic theoretical description of the relation between V" and doping exist. 
For example the V" formation energy, in presence of doping, is usually considered™ only 
for the V" to dopant atoms ratio, y v ** / x , equal to 1 and again V" are assumed to induce 
delocalized states which could mediate the electron conduction in case of resistive switching, 
regardless of the value of y v ^ j x . 

Instead, in case of X +3 elements, like iron, the most stable configuration is expected 
to have y v ^ / x = 0.5 for charge compensation^. We thus focus our attention on this 
configuration describing how the properties of the system would change if y v ** j x deviates 
from the value 0.5. 

In sec. Uwe describe both the theoretical and the experimental approach to the description 
of ZrO^-Fe. The results from first-principles simulations are presented in sec. [Ill The 
electronic and structural properties of the system are described in function of the doping 
and oxygen vacancies concentration within density functional theory (DFT) in the standard 
generalized gradient approximation (GGA). For TM oxides the standard approximations to 
DFT are known to fail in the description of the so called on-site correlation. Thus DFT can 
be corrected with a "Hubbard" term, DFT+[7 scheme, where U is an external parameter, 
which improves the DOS of the valence electrons. However little is known in the case of 
TM used as dopant in DMO. Thus we also investigate how this term would influence the 
electronic properties of the system in ZrO^-Fe. 

The experimental results are then presented in sec. IIHI Here we show that, indeed, the 
measured properties best agrees with the y v ** j x = 0.5 configuration. Moreover a detailed 
comparison of the measured valence band (VB) and DFT density of states (DOS) is done. 
This is a direct way to explore the value of the on-site electronic correlation on this system, 
i.e. to adjust the value of the U parameter to be used in the DFT+?7 approch. 
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I. FRAMEWORK 



A. Computational approach 

We computed, from first-principles, the ground state of the two most common phases 
of ZrOii i- e - t ne tetragonal and the monoclinic phases, at different doping concentrations. 
We used the PWscf (4.3.2) package^, considering a super-cell with 96 atoms (few less 
when V" are considered) and in some cases also a smaller super-cell with 12 atoms for the 
description of the highest doping configuration. For all systems the atomic positions are 
fully relaxed. The ground state was computed within the GGA— to the DFT scheme^^ 
with ultra-soft pseudo-potentials^ 1 ^. We used a 35 Ry cut — off for the wave-functions, 
400 Ry cut — off for the augmentation density and a Monkhorst-Pack grid 2x2x2 for the 
Brillouin zone to have the error on the energy differences between the monoclinic and the 
tetragonal phase lower than 1 meV per formula unit (f.u.); this was the most stringent 
condition for our simulations. We estimated the error on the total energy to be lower 
than 0.1 eV / f.u.. Convergence paramenters are 10~ 8 Ry on the total energy for the scf 
cycles and both 10~ 4 Ry on the total energy and 10~ 3 Ry/Bohr on the forces for the 
atomic relaxation. The pseudo-potential of Zr includes semi-core electrons. Fe atoms 
were placed at the substitutional Zr sites and kept as far as possible from each other to 
mimic uniform doping. For V" instead we considered many different configurations (see 
discussion in sec. [TTJ) , specifically we considered ZrO^-Fe at the atomic doping concentration 
xp e = 6.25%, 12.5%, 18.75%, 25% with, y V " /Fe = an d without, y V " /Fe = 0' oxygen 
vacancies. We also considered y v " /Fe = 1-0 f° r x Fe = 6.25%, 25%. In total we studied 
about 50 different systems of Zri^ x Fe x 02~ z V" z changing %Fe and z v »* = ip e x y v »* / Fe 
for either the monoclinic or the tetragonal structure. For few selected configurations, i.e. 
at the lowest and the highest considered doping concentrations xp e = 6.25, 25%, we also 
performed calculations within the simplified GGA+[7 approach 44 implemented in the PWscf 
package, again considering y v ** j Fe = 0, 0.5, 1, in order to explore the effect of the Hubbard 
correction on the electronic structure of the system. The results are presented mainly for the 
high-doping situation which we have also experimentally. The configurations at y v »* j Fe = 
and 1 resulted to be metallic and in these case the convergence of the physical quantities 
against the sampling of the k-points grid was verified. 
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The cell parameters for both the tetragonal and the monoclinic phase of pure ZrOi are 



the same used in Ref. 



31 



Specifically for the monoclinic phase a = 5.18 A, b/a = 1.011, 
c/a = 1.037 and /3 = 99°10'; while for the tetragonal phase a = 5.18 Aand c/a = 1.0305. 
The same parameters were used for ZrO^-Fe as well. However we even performed a full 
relaxation of our 96 atoms super-cell for few selected configurations and we found out that 
this have a negligible impact on the properties of the system here considered. 

In sec. [TTl we systematically compare the results of the present simulations with the 
ZrO^y (Y doped ZrO?) system. Yttrium is one of the most studied and used dopant 
of ZrC>2 and shares with iron the same valence. All the data reported for ZrO^-Y are from 
ref. 311 



In order to describe the semi-core levels of iron and compare the results with XPS mea- 
surements, we run calculations with a norm-conserving fully-relativistic approach. To this 
end, we used Hartwigsen, Goedecker, and Hutte (HGH) pseudo-potentials^, which con- 
tain semi-core electrons in valence and are constructed with a fully relativistic calcula- 
tion. The latter are not available within the PWscf 2 - code and so we used the abinit (6.8) 
code~ We studied the semi-core levels only for the xp e = 25% at. case again considering 
Vv" /Fe = 0) 0-5; 1- We used smaller super-cells, 12 atoms i]Jv" /Fe = an d 1) an d a 24 
atoms supercell (Vv" /Fe = 0-5), with cut-off of 170 Ry and a Monkhorst-Pack grid 3x3x3 
and 3x3x2 respectively for the Brillouin zone to have the error on the energy levels position 
lower then 0.1 eV. The very high energy cut-off was needed, as the norm conserving HGH 
pseudo-potentials are harder than the ultra-soft ones used with PWscf and also because the 
semi-core levels are much more localized than valence electrons. The value xp e = 25% at. 
was chosen to have smaller super-cells but also because this is quite close to the experi- 
mentally measured doping concentration in our films. The atomic positions instead were 
obtained relaxing the same structures with the PWscf code and then we checked that the 
residual forces on the atoms computed with Abinit were negligible. 

Finally for a quantitative comparison of the measured photo-emission and the computed 
valence band we have performed calculations within GGA+[7 at U = 1.0, 2.0, 3.3 eV at 
%Fe = 18.75% at. and y v ** j Fe = 0.5. A theoretical smearing of 0.02 Ry was used to generate 
the DOS used in Figs. [SHwhile a higher smearing of 0.06 Ry was used for the DOS in Fig. [7] 
to mimic the experimental peak width. 



B. Experimental setup 



Experimentally Zr0 2 and Zr0 2 :Fe thin films were grown on Si/Si02 substrates in a 
flow-type hot wall atomic layer deposition reactor (ASM F120) starting from /3-diketonates 
metalorganic precursors, namely Zr{C\\H 1902)4 for Zr and Fe(CnHi 9 2 ) 3 for Fe. To grant 
a stable reactivity, Zr precursor was kept at 170°C, while Fe precursor was maintained at 
115°C. Ozone was used as oxidizing gas in the reaction process The film growth was achieved 
by alternately introducing the reactants separated by N 2 inert gas purging pulses. The Fe 
concentration in Zr0 2 :Fe films was tuned tailoring the Zr/Fe precursors pulsing ratio and 
the growth temperature was maintained at 350°C (details in Ref. I29L After the deposition 
the films were annealed at 600° C in N 2 flux for 60s. The growth parameters were tuned in 
order to fix the thickness, d — 19 ± 1 nm, and the doping concentration xp e = 20% ± 3% 
for the Zr0 2 :Fe films. XFe was chosen in order to stabilize the tetragonal phase according 
to our theoretical results. 



Film crystallinity was checked by X-ray diffraction (XRD) at fixed grazing incidence 
angle oj — 1° and using Cu K a (A = 0.154 nm) monochromated and collimated X-ray 
beam (details in Ref. l30l ). Film uniform doping along its thickness was checked by Time 
of Flight Secondary Ion Mass Spectrometry (ToF-SIMS) depth profiling using an ION 
TOF IV instrument, with 500 eV Cs + ions for sputtering and 25 keV Ga + ions for analysis. 
Secondary ions were collected in negative polarity and interlaced mode. Recorded intensities 
were normalized to 30 Si intensity in bulk silicon. The instrument depth resolution is below 
1 nanometer. 



To elucidate Fe chemical state and concentration in Zr0 2 :Fe films, X-ray photo- 
emission (XPS) measurements were performed on a PHI 5600 instrument equipped with a 
monochromatic Al K a x-ray source (E = 1486.6 eV) and a concentric hemispherical ana- 
lyzer. The spectra were collected at a take-off angle of 45° and band-pass energy 11.50 eV. 
The instrument resolution is 0.5 eV. 
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II. FIRST PRINCIPLES PREDICTIONS 



A. V" and structural properties 

In the literature Zr0 2 :Fe has been studied as a candidate material for oxygen sensing 
applications because Fe +3 atoms, replacing Zr +A atoms, are expected to induce oxygen 
vacancies for charge compensation^*™ Thus, as a first step, we consider the V" formation 
energy: 

AE 1 (x Fe , z) = ( E[Zn. x Fe x 2 } - ( E[Z n _ x Fe x 2 . z } + {z/2) t i[0 2 \ ) )/z (1) 

at fixed y V " /Fe — z/x — 0.5, i.e. for a charge compensated system. Here we considered 
both the oxygen rich condition (Fig. HJ(a), n[0 2 ] = -^[Cy with E^Cy the total energy of 
an isolated oxygen molecule in its ground state) and the oxygen poor condition (Fig. [TJ (b) , 
^[0 2 ] = E[Zr0 2 ] — E[Zr]). The formation energy for Zr0 2 :Fe is compared with the case 
of pure Zr0 2 , AEi(0,z) and Zr0 2 :Y, AEi(xy, z). To this end we considered different V" 
concentrations and, for each, different V" configurations. However we found that AE(x, z) 
is mainly determined by the kind of dopant, while the influence of the other parameters 
is lower. In Fig. [T] the changes due of these parameters results in different values for each 
system. While the Vq formation energy is negative in Zr0 2 :Y already in the oxygen rich 
case, in Zr0 2 :Fe films it is slightly positive, i.e. AE{ etra m 0.5 eV, but ten times lower 
than in pure Zr0 2 . Varying the chemical potential from the oxygen rich to the oxygen poor 
configuration AE\ etra becomes negative, thus Fe favors the formation of V" . 

The creation of oxygen vacancies induces disorder in the system (see also the inset in 
Fig. ED thus the most symmetric phases are expected to be favored against the monoclinic 
phase. To evaluate this effect in Fig. [2] we consider the energy difference AE 2 (x) between 
the tetragonal and the monoclinic phase as a function of the doping concentration, at fixed 
Vv" jx = 0-5. We look for the iron atomic percent, Xp e , at which the tetragonal phase 
becomes favored. 

The value of AE 2 is very small and thus at the limit of the DFT-GGA resolution. The 
computed energy difference between the two phases at zero doping is AE 2 (0) = 109 meV/ f.u., 
in agreement with previous works, (63 meV/f.u.—, 144 meV/ /.«.—); the experimental esti- 
mation is 63 meV/f.u.—. It is reasonable to assume that the trend of the energy difference is 
better computed than its absolute value and accordingly, assuming a constant "zero-doping 
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oxygen rich oxygen poor 




Figure 1. (color online) DFT (GGA) formation energy of oxygen vacancies (see Eq. [T]) in Zr0 2 '.Y, 
Zr0 2 :Fe and Zr0 2 in the two extrema case of (a) oxygen rich conditions and (b) oxygen poor 
conditions in both the tetragonal and the monoclinic structure. The doped systems are considered 
in the charge compensated configuration (i.e. Uv"/x = ^-5 for X = Fe,Y). The values are 
computed for different oxygen vacancies concentrations and also varying, for some concentrations, 
the position of the oxygen vacancies. In panels (a) and (b) histograms are presented in the same 
order (and colors). 

error" of ~ 46 meV/f.u. for every Fe concentration, we can subtract it. Being AE 2 of the 
order of few meV/ f.u. also the phonon energy of the two lattice could play a role. Indeed the 
monoclinic to tetragonal phase transition at ~ 1440 K can be explained in this terms^^i. 
Thus we considered the energy difference of the lattice between the two structures for the 
undoped system. At room temperature however we found this contribution to be almost 
negligible, ~ 5 meV/ f.u.. 

AE 2 come out, instead, to be particularly sensitive to the chosen atomic configuration. 
Accordingly the data in Fig. [2] are scattered, with AE 2 changing of few meV/ f.u. at given 
xp e . To extract the exact Xp e a statistical occupation of the different configurations should 
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Figure 2. (color online) DFT (GGA) total energy difference per formula unit between the tetragonal 
against the monoclinic phase for ZrO^-Fe, panel (a), and ZrOi'Y , panel (b). Total energies are 
computed for the charge compensated system (dots, Uy" /x = ^-5 for X = Fe,Y) changing the 
atomic configurations for each given concentrations. The shadowed areas are guides for the eyes 
while the continuous lines are a linear fit of the data. Also the results for the systems without 
oxygen vacancies (crosses, y v " j x = 0) are shown for comparison. The zero level is shifted of 
(i) —46 meV/ f.u. to align the energy difference at zero doping with the experimental value, (ii) 
—5 meV/ f.u. to include the computed zero-point-energy difference of the two lattices. 

be considered. However, to this end, one should sample a huge number of configurations, 
which is not feasible within DFT. In the present paper we assumed that, fixed xp e and z v >» , 
changing the configurations for the V" , AE 2 spans uniformly a given energy range (Er) 
which can be extrapolated considering a limited number of configurations. Er is expected 
to increase, increasing the doping concentration, as an increasing number of configurations 
becomes available. With this assumptions x% e was extracted considering the central value 
of the Er. 

In practice this was done with a linear fit of the data. In Fig. [2J to obtain the critical 
doping concentration, AE 2 at zero doping is matched at the experimental value 63 meV, 
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while the theoretical results would be 109 meV. The result, x Fe ~ 12% at., can be compared 
with the case of ZrO^'Y ', where the same approach gives x Y ' ~ 7% at. which exactly matches 
the experimental valued. We stress that with this approach the exact doping concentration 
can be affected by an error which can be as large as few atomic percent. What is significant 
here is the comparison of the two systems, i.e. ZrO^Y and ZrO^.Fe. Indeed both dopants, 
inducing oxygen vacancies favor the tetragonal against the monoclinic structure. However 
the two linear fits posses different slopes and we can conclude that iron is less efficient than 
yttrium in inducing a monoclinic to tetragonal phase transition. In Fig. [2] we also report 
the energy difference between the monoclinic and the tetragonal phase for the case without 
oxygen vacancies, i.e. Vy" /x = 0- ^ n this configuration we found that the local structure 
of the crystal is much less distorted by doping and accordingly the variation of the energy 
difference between the two phases is small. This confirms that a key role in the monoclinic 
to tetragonal phase transition is played by oxygen vacancies^ and not by the dopant itself. 

B. Electronic properties 

Given the results of the previous section and the fact that experimentally we describe a 
system at high doping concentration, which we found to be in the tetragonal phase, in the 
description of the electonic properties of the system we focus our attention on the tetragonal 
structure of ZrO^-Fe. 

The main difference between Y and Fe is the presence of the unfilled Fe(d) orbitals which, 
falling inside the energy gap of zirconia, determine the electronic properties of the doped 
system. The d-orbitals occupation is also strongly affected by V" and is used here, together 
with the computed magnetic moment, to infer the Fe oxidation state. At y V '» / Fe = 0, Fe 
acts as an acceptor (Fig. 0(a)) with the creation of holes in the majority spin VB. These 
are preferentially located on the Fe(d) orbitals as shown by the projected-DOS with the 
projection of the hole states on the d-orbitals close to 0.5. Thus Fe is forced in the Fe +A 
oxidation state with a magnetic moment per atom equal to 4 Bohr magnetons (/is). The 
creation of V" release the electrons captured by the O anions. At y v ** ^ Fe = 0.5 the system 
turns into a charge-transfer semi-conductor (see Fig. [31(b)), i.e. the V" do not create an 
impurity bands, as it would happen in Zr0 2 , but compensate the holes in the Fe(d) orbitals. 
In this configuration Fe atoms are in the +3 oxidation state and the magnetic moment per 
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Figure 3. (color online) Total (full line) and d-orbital projected (dashed line) density of states 
(DOS) at the GGA level of ZrO^'-Fe. at XFe = 25% with y v »» ^ Fe equal to respectively (panel 
a), 0.5 (panel b), 1 (panel c). The vertical dashed line marks the Fermi level. The Fermi level of 
panel (b) is the zero of the energy axis, while in panels (a) and (c) the zero is obtained aligning 
the bottom of the valence band at ~ —6.5 eV as in panel (b). 

iron atom is maximized, 5 fis- If Vv" /Fe exceeds 0.5, electrons start to fill the minority 
Fe(d) levels. This decreases the average magnetic moment, while the system reverts to an 
half-metal. At y v ** j Fe = 1 (Fig. EJ(c)) all iron atoms are in a +2 oxidation state with the 
per atom magnetic moment equal to 4 In Fig. [3]we also notice that at y v ** j Fe < 0.5 no 
extra state, other than the Fe(d) orbitals, appears between the valence and the conduction 
band of Zr0 2 - Only when y v »* j Fe > 0.5 (Fig. EJ(c)) such a state exists. The latter can be 
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associated to an impurity band which has been suggested to create bound magnetic polarons 
in case of magnetic doping-. However the configuration y v ^ ^ Fe = 1 is not favored. Indeed 
the energy cost, for each extra V" created in the system, of the reaction 

Z ri _ x 2 _ x/2 Fe x -+ Zr^ x 2 „ x Fe x + ^[0 2 ] (2) 

changes from w 2.5 eV, oxygen rich conditions, to w 0. eV, oxygen poor conditions, thus 
remaining positive for any value of the oxygen chemical potential. As for the V" formation 
energy AEi, this value is weakly dependent on the atomic doping xf s . Last but not least, 
even if in this case the impurity band exist, it is empty. Thus the possible existence of 
bound magnetic polaron in Zr0 2 :Fe is unlike. We will also show in the next section that, 
experimentally, iron in Zr0 2 :Fe is in the Fe +3 and not in the Fe +2 oxidation state. 

We remark that, even if at y v »» / Fe = and y v ** j Fe = 1 the system is metallic, the per 
atom magnetic moment is integer. The reason is that in both cases Zr0 2 :Fe is indeed an 
half-metal and thus electrons can move across the Fermi level only in one spin channel. We 
have verified this result increasing the sampling of the k-points grid from 2x2x2 to 3x3x3 
in the 96 super-cell and from 4x4x4 to 8x8x8 in the 12 atoms super-cell. In both cases 
the system remains metallic, with fractional occupation in the majority {Uv" /Fe = 0) or 
minority (Vy" /Fe = 1) s P m channel (a smearing of 0.002 Ry was used in the self-consisten 
cycle in this case), but with constant per atom magnetic moment m z = 4 /xg. In principle 
the 96 atoms super-cell with a sampling 3x3x3 is equivalent to the 12 atoms supercell with 
sampling 6x6x6. However the two could differ because in the 96 atoms super-cell, removing 
symmetries, disorder is taken into account. This could for example induce a localization of 
holes on the Fe atoms. Thus the convergence check were also a rough way to explore possible 
Anderson-like localization mechanisms. However we did not observe such phenomena. 

These are the prediction of the GGA. Howerver for TM oxides this approximation is 
known to suffer of some deficiencies. In particular it suffers of the well known problem of 
self-interaction, which tends to delocalize too much the d orbitals. A common way to avoid 
this problem is to correct the DFT scheme with a Hubbard like term U which enters as an 
external parameter. The value of U is system dependent and should be optimized either 
with a direct comparison with experimental data or with a self-consistent approach. In the 
literature usually U — 1 — 3 eV for elemental iron and [7 = 2 — 6 eV in iron oxides. For 
example Cococcini et al.— report, after a self-consistent calculation, U ~ 2.2 eV for metal 
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Figure 4. (color online) Total (full line) and d-orbital projected (dashed line) density of states 
(DOS) at the GGA+[7 level of ZrO^-Fe at xj? e = 25% with y v »» j Fe equal to respectively (panel 
a), 0.5 (panel b), 1 (panel c). The vertical dashed line marks the Fermi level. The Fermi level of 
panel (b) is the zero of the energy axis, while in panels (a) and (c) the zero is obtained aligning 
the top of the conduction band at ~ 5 eV as in panel (b). 

iron and U = 4.3 eV for FeO. Here we begin choosing an intermediate value, U = 3.3 eV, 
in order to evaluate the physical effects introduced by this correction. 

In Fig. S] we plot the DOS for the GGA+U approach at y v " /Fe = 0' 1- We can 
directly compare the results with the GGA DOS plots in Fig. [31 As expected the U correction 
pushed down the occupied d level and a sharp structure appeared in the DOS just below 
the VB of ZrO-i.Fe. Also the crystal field splitting of the spin minority d orbitals, between 
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the e g and the t g states, was reduced, and is not distinguishable anymore with the smearing 
parameter used in the plot; with the exception of the case y v »» j Fe = 1. However in the 
charge-compensated situation, y v « j Fe = 0.5, these corrections do not alter the qualitative 
description of the system, which remains a magnetic semi-conductor with the magnetic 
moment per atom maximized. Instead, when we deviate from this configuration, we notice 
two important differences. For y v '» j Fe < 0.5 the holes created in the VB are less localized 
on the Fe atoms. Indeed the projection of the hole states on the d-levels drops from « 0.5 
(GGA) to less than 0.1 (GGA+?7). Thus iron is in the Fe 3+ configuration, while the holes 
are in the ZrC>2 VB, i.e. on the oxygen atoms. Accordingly the V" formation energy drops 
from 0.5 eV (GGA) to 0.0 eV because oxygen atoms are more weakly bound to the system. 
For y v »* j Fe > 0.5 the extra electrons start to fill the minority (i-levels, as in the GGA case. 
However the newly occupied levels are pushed down in energy and thus the system is not 
metallic but it displays an energy gap, i.e. GGA+[7 predicts a a Mott insulator in this case. 
Also for the GGA + U case we verified that in the metallic case (i.e. at y v ** j Fe = 0) the 
value of the magnetic moment remains constant improving the sampling of the Brillouin 
zone. 

The electronic properties in the present section were reported for xp e = 25%. We did not 
find significant changes for the other doping concentrations, at least for y v ** j Fe = 0.5. At 
the lowest computed doping concentration however, xp e = 6.25%, the Fe atoms are too far 
apart and the localized rf-levels do not create a band. Thus the metallic phases predicted 
within GGA (y v " /Fe = 0> -0 become semi-conducting phases with defects states localized 
close to the Fermi level. 

III. EXPERIMENTAL RESULTS 
A. Structural characterization 

Experimentally, as a first step, we first studied the structural properties of the films 
growth by atomic layer deposition. In Fig. [5] the ToF-SIMS depth profile of a representa- 
tive film (namely, ZrO^-Fe at Xp e — 20% at.), including Fe, FeO, ZrO and Si negative 
secondary ion intensity profile is graphed. Fe and FeO are both used as representative 
of Fe distribution along the film thickness; in particular FeO ion fragment has not to be 
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Figure 5. (color online) Tof-SIMS depth profile of Zr0 2 :Fe at x Fe « 20% at. . 

considered as a mark of FeO chemical compound in the film, but as a fingerprint of Fe 
embedded in the Zr0 2 host matrix. The flatness of ZrO and Fe related profiles indicates 
that the film grows uniformly during the ALD process, without changes in the distribution 
of the chemical species, evidencing that the growth process is well controlled. Further, Si 
diffusion in Zr0 2 is excluded with a well distinct film/substrate interface, an indication 
that the substrate does not affect Zr0 2 :Fe properties both during the film growth and the 
thermal treatment. Furthermore, the Fe profile is almost constant, thus it is the doping 
in the film, and the absence of large fluctuations such as peaked maxima, can exclude Fe 
clustering. Indeed the latter would have been observed as a sudden increase of Fe intensity 
with a concomitant abrupt decrease of FeO intensity, indicating that an Fe rich / O poor 
environment is detected. Instead both Fe and FeO signals mimic the same profile shape, 
confirming that Fe is uniformly diluted within the Zr0 2 matrix. 
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Figure 6. (color online) - XRD patterns of ZrO<i (blue) and ZrOi'-Fe (red) (Fe doping ~ 20% at.) 
films evidencing Fe doping is effective in suppressing the monoclinic phase. {.{ZrO-i) and m(Zr02) 
indicates the reflections from reference tetragonal and monoclinic Zr02, respectively 3 ^. On the 
right the relaxed DFT structure for m(ZrOo)) XiZrO-i) and \,{ZrOi'-F e) at XFe = 25% at. repre- 
sented with the xcrysden package (see Ref . I39I) ; Zr atoms in blue, Fe atoms in red and the smaller 
O atoms in black. 

To get details on the film crystalline structure, in Fig. [6] we compare the XRD patterns of 
ZrO<i and ZrO^-Fe. Both films mainly present the cubic/tetragonal phase. Indeed in these 
films there is a balance between the bulk energy, where the monoclinic phase is favored, 
and the surface energy, where the tetragonal phase is favored. The critical grain size 3 ^* 3 ^ 
below which the tetragonal phase become the most favored is ~ 15 nm. In our films, being 
the grain size close to the film thickness (from XRD data), we are close to this critical 
value. This can be evinced from the XRD patterns of pure ZrOi where the peaks of the 
monoclinic phase are also evident. However in the ZrO^'Fe films the monoclinic phase is 
completely suppressed, confirming our theoretical findings. Even from these measures there 
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Figure 7. (color online) (a): The Fe(2p) core level photo-emission spectra in ZrC>2'-Fe. (b): 
Fe(3p) and Zr(As) photo-emission spectra and computed DOS for ZrOi'-Fe with y v '» ^ Fe equal to 
(green dashed), 0.5 (red continuous), 1 (black dot-dashed). The Fe(3p) majority spin level is in 
light gray, (c): Measured and computed valence band (VB) for pure ZrOi- (d): Measured VB for 
ZrO<i'-Fe. Computed VB for ZrOi'-Fe with Fe doping substitutional at y v " j Fe = 0.5 (continuous 
red line) or interstitials (dashed maroon line). In panels (b)-(c)-(d) the experimental data (and 
fit) are vertically shifted respect to the DFT-DOS. All DOS are obtained at within the GGA. The 
experimental data were collected with the PHI 5600 instrument (see details in sec. II Bp 

is no indication of segregated iron phase or iron oxide clusters. 
B. Electronic properties 

In Fig. [7] we report the high resolution spectra of the Fe(2p) core level (a), the Fe(3p) 
semi-core (b) levels and the VB (c-d). In Fig. [7|(a-b) the data were fitted with a doublet of 
asymmetric Voigt functions for the two main peaks plus a Voigt function for the satellite on 
top of a Shirley background and in Fig. El(b-d) the spectra are compared with DFT(GGA)- 
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Table I. Energy distances [eV] from Fe(2p 3 / 2 ). 
Data for iron oxides from Ref. 43 



Fe 2 O s Fe 3 0± FeO Zr0 2 :Fe 



Fe(2p 1/2 ) 


-13.6 


-13.5 


-13.6 


-13.5 


satellite 


-7.8 


not pres. 


-6.0 


-8.6 



Fe(3p) 655.4 not av. 653.9 655.2 



DOS computed as described in sec. II Al 

The change of the XPS-VB from Zr0 2 (blue) to Zr0 2 :Fe (red) is in agreement with 
the DFT(GGA)-DOS obtained considering substitutional iron doping. In particular exper- 
imentally the double peak structure of pure Zr0 2 is suppressed with doping. Theoretically 
this behavior is reproduced only assuming substitional doping. 

The core or semi-core levels of TM usually show a structured shape due to, at least, 
four factors: the spin-orbit (SO) splitting, the exchange splitting, the multiplet splitting 
and the eh screening to the core-hole. The SO term is responsible for the 2pi/ 2 - 2p 3 / 2 
splitting AEso — 13.5 eV and is not sensitive to the chemical environment (see Fig. 13(a)). 
The exchange and multiplet splitting instead give the characteristic asymmetric shape of the 
XPS peaks in metals. Finally the screening effect, which is strongly sensitive to the chemical 
environment^"—, can create satellites. For the Fe(2p) core level the distance between the 
satellite and the Fe{2p 3 / 2 ) peak is a marker of the iron oxidation stated. Also the position of 



the Fe(3p) peak (Fig. 0fb 



43 



I is sensitive to the Fe chemical environment^. The comparison 
reported in Table [U shows that iron is in the Fe +3 oxidation 



with the values of Ref. 
state. 

According to our DFT results the Fe oxidation state is strongly related to the presence 
of V in the system (see Fig. [3]). To better describe this point we study the Fe{3p) semi- 
core levels with first principles simulations. Indeed the Fe(3p) wave-functions are spatially 
localized close to the Fe(3d), which are in valence, and so are very sensitive to the chemical 
environment. The energy of the Zr(4s) level is used as a reference to properly align the 
experimental XPS levels with the theoretical DOS. 

In our approach the SO coupling term was included, both in the pseudo-potentials and 
in the hamiltonian, while the multiplet and the exchange splitting were accounted for by 
the exchange-correlation (xc) potential. For the Fe(3p) level we found AEso < 1 eV, while 
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AE XC 5 eV between the spin minority and the spin majority which is clearly visible in 
Fig. 0(b). This is overestimated by DFT. In the case of semi-core levels Takahashi et al. 41 
showed that the screening effects, which are not included in the present approach, gives a 
broadening and a shift of the majority spin channel with, possibly, the creation of satellites. 
Indeed we can suppose that these effects would correct the overestimated AE XC fh 5 eV, 
shifting the majority-spin energy level close to the minority one giving a single asymmetric 
peak with higher intensity as in the experimental case. However such an approach is beyond 
the scope of the present work. The minority-spin channel instead is less affected by screening 
effects retaining the independent-particle structure with the onset of the spectrum due to 
absorption from this channel 41 -. Thus we compared the energy position of the minority 
DOS with the measured Fe(3p) XPS spectrum. In our simulations the distance of the 
Fe(3p) minority peak from the Zr(4s) level, AE y , is strongly dependent on y v »* j Fe with 
AE y=Q = 1.8, AE y= o. 5 = 3.1 and AE y =i = 4.0 eV. The value AE y=Q5 , i.e. the configuration 
with iron in the Fe +3 oxidation state, best agrees with the experimentally measured splitting 
AE = 2.9 eV, in agreement with the conclusion drawn from Table [I] and in general from 
seaHl 

In sec. [TT] we showed that the electronic properties and in particular the shape of the 
valence band could be strongly influenced by the on-site electronic correlation, by the com- 
parison of GGA and GGA+U predictions at U = 3.3 eV. 

To decide which of the two scenarios, GGA or GGA+U, better describes the experimental 
situation we compared the obtained DOS for the charge-compensated case with the mea- 
sured photo-emission from the VB; the values U = 0.0, 1.0, 2.0, 3.3 eV are here considered. 
To mimic the experimental spectrum, we have superimposed to the DFT-DOS a Shirley like 
background, i.e. a background proportional to the integral of the DOS. Also, for a quanti- 
tative comparison, we computed the theoretical DOS for y v »* j Fe = 0.5 and xp e = 18.75%, 
which is the theoretical value closest to the experimental measured doping. 

In Fig. [8] we see that the structure which identifies the d levels in the GGA+{7, at the 
reference value U = 3.3 eV, is not present experimentally and the agreement between theory 
and experiment is much better in the standard GGA (i.e. U = 0. eV). At the intermediate 
values U = 1.0 eV and U = 2.0 eV such structure is not visible, however the agreement 
with the experimental results is worse than for the {7 = 0. eV case. We can conclude that 
the value {7 = best agrees with the photo-emission VB, and that, given the experimental 
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Figure 8. (color online) Valence band of Zr02-Fe. The GGA+U scheme at xpe = 18.75 %, 
y v »» j Fe = 0.5, for the values of U = 0.0, 1.0, 2.0, 3.3, is compared against experimental data. 
The smearing parameter used for the plot is 0.06 Ry. The experimental data were collected with 
the PHI 5600 instrument (see details in sec. II Bp . 

resolution, the optimal choice of U must be between and 1 eV. Thus in ZrO^-Fe the 
effect of the self-interaction of the d orbitals, which is corrected by the Hubbard U term, 
is smaller than in common iron oxides. This is an "a posteriori" justification of the results 
obtained, in the present work, within the GGA. 



IV. CONCLUSIONS 

In conclusion we studied iron doped zirconia both theoretically, with first-principles sim- 
ulations, and experimentally, with structural, chemical and electronic characterization of 
thin films grown by atomic layer deposition. 

As expected from simple considerations, iron was found experimentally in the Fe +3 ox- 
idation state. We also found that it induces a monoclinic to tetragonal phase transition. 
Theoretically the oxidation state was related to presence of oxygen vacancies which play a 
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key role in the structural phase transition. The theoretical findings have been tested with 
a detailed comparison against photo-emission spectra of the samples grown by atomic layer 
deposition to validate the assumptions. These results are a confirmation that iron doped 
zirconia could be a good candidate in view of oxygen sensing applications as reported in the 
past. 

Moreover the presence of vacancies is seen not only to influence the structure of the 
system but, theoretically, also to determine the density of states at the Fermi level and 
the eventual presence of impurity states in the gap which could be associated to magnetic 
polarons. In particular, we discussed how the ratio between oxygen vacancies and the 
iron atoms concentration shifts the Fermi level of the system. We found that in the most 
stable configuration, the Fe +3 iron atoms are charge-compensated by the presence of oxygen 
vacancies with a ration of 0.5, i.e. one vacancy each two iron atoms. The resulting system 
is a semi-conductor with no impurity state in the gap. 

These results should be considered for a correct description of the behavior of iron doped 
zirconia, or more in general of high-k oxides doped with valence +3 elements, in resistive 
switching devices. Moreover the absence of impurity states rules out the magnetic polaron 
model as a possible mechanism to explain the magnetic properties of the system. 

Finally we have explored the importance of the Hubbard U correction. Indeed, theoret- 
ically, varying the value of U from eV to 3.3 eV the electronic propertes of the system 
change significantly. We showed that in iron doped zirconia the value U w eV best agrees 
with the experimental data, thus indicating that the on site electronic correlation is low in 
this system. 
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